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ABSTRACT 

We present new results from the Disks@EVLA program for two young stars: CY Tau and DoAr 25. We 
trace continuum emission arising from their circusmtellar disks from spatially resolved observations, 
down to tens of AU scales, at A = 0.9, 2.8,8.0, 9.8 mm for DoAr 25 and at A = 1.3, 2.8, 7.1 mm for 
GY Tau. Additionally, we constrain the amount of emission whose origin is different from thermal 
dust emission from 5 cm observations. Directly from interferometric data, we find that observations 
at 7 mm and 1 cm trace emission from a compact disk while millimeter-wave observations trace an 
extended disk structure. From a physical disk model, where we characterize the disk structure of 
GY Tau and DoAr 25 at wavelengths shorter than 5 cm, we find that (1) dust continuum emission is 
optically thin at the observed wavelengths and over the spatial scales studied, (2) a constant value of 
the dust opacity is not warranted by our observations, and (3) a high-significance radial gradient of 
the dust opacity spectral index, /S, is consistent with the observed dust emission in both disks, with 
low-/3 in the inner disk and high-/? in the outer disk. Assuming that changes in dust properties arise 
solely due to changes in the maximum particle size (ttmax), we constrain radial variations of Umax in 
both disks, from cm-sized particles in the inner disk (R < 40 AU) to millimeter sizes in the outer disk 
(i? > 80 AU). These observational constraints agree with theoretical predictions of the radial-drift 
barrier, however, fragmentation of dust grains could explain our amax[R) constraints if these disks 
have lower turbulence and/or if dust can survive high-velocity collisions. 

Subject headings: stars: individual (GY Tau, DoAr 25) — protoplanetary disks — stars: formation 
— radio continuum: planetary systems 


1. INTRODUGTION 

The process of planet formation requires that small 
dust grains in the interstellar medium (ISM) undergo a 
dramatic transformation, growing by more than 12 or¬ 
ders of magnitude in order to form terrestrial planets 
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and the cores of gas and ice giants. As a very first step, 
these /im-sized ISM dust g rains must increas e their size 
and become macroscopic (Testi et al. 20141, a process 
that alters the optical properties of the dust grains con¬ 
siderably: as grains reach millimeter or larger sizes, the 
absolute value of the dust opacity, kx, decreases, and 
at the same time the power-law spectral index of the 
dust opacity, S (where k\ oc A~^), becomes smaller (e.g. 
Miyake fc Nakagawa||1993[ [Henning fc Stognienko||T99^ 
Draine||2006| . Since thermal dust emission at millimeter 
and centimeter wavelengths is (mostly) in the optically 
thin regime, the observed dust continuum emission will 
trace the bulk of the disk mass modulated by the dust 
opacity: S,^ oe kxxTjX B^{T), where S is the dust mass 
surface density and B^{T) is the Planck function evalu¬ 
ated at the dust temperature T. For optically thin dust, 
warm enough to be in the Rayleigh-Jeans (R-J) regime 
{hv kT), it follows that oc Thus, direct 

measurements of the dust emission spectrum (i.e., the 
spectral energy distribution, SED, at long wavelengths) 
can be used to derive the value o f /?, a method exten¬ 
sively employed in t he literature ( Beckwith&_Sarge^ 


1991 Henning et al. [19^ [Testi et al.|20UH [Galvet et al. 


Testi et al.[[20d3[[Natta & Testi|2004| [ Wilner et al. 

” ^ 'J Will . 


2002 


2005 Andrews & Williams! 


odmann et al.[[2006r 


Andrews fc Williams|2007| [Lommen et al.[[2009[ [Ricci et| 


al.[[2010a|b[[20111 [U bach et al.[[2012D . 

Disk-integrated measurements, from sub-mm to cm 
wavelengths, have shown that in most protoplanetary 
disks the value of f3 tends to be lower than in the ISM 
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(where /3ism 1-5 — 2.0, consistent with the presence 
of /rm-siz ed dust grains or smaller in the interstellar 
medium, Li & Draine||2001 ). 

Possible explanations for the low values of /3 include 
the presence of regions of high optical depth, which will 
drive the dust emission spectrum to be close to oc A“^, 
and t hus a low 0 would be inferred (incorrectly). How¬ 


ever, 


Ricci et al. 


(] (2012|) has shown this would be only 
"brightest and most massive disks. Dif- 


plauslble tor the 
ferent composition and/or porosit y could affect the opti¬ 
cal properties of dust grains (e.g., Henning fc Mutschke 
19971 [Semenov et al.||2003[ Houdet et al.||2UU5| |Katao|a 
et al. |2014p , but tor a sensible set of dust properties the 
most signiiicant influence on the o pacity spectra l index 
arises from an increase in grain size (Draine|2006 ). Thus, 


the difference in the dust opacity spectral index of disks, 
as compared to the ISM, implies that growth of at least 4 
orders of magnitude in size has t aken place inside circum- 
stellar disk s (for a review see: Natta et al. 2007 Testi 
et al. 2014). But since most of these measurements are 
of the global disk properties, they provide limited knowl¬ 
edge of the evolution of the dust properties as a function 
of the distance to the central star (hereafter referred to 
as radius). 

Theoretical analyses predict that the average dust 
grain size will change with radius, as dust grows an d frag¬ 
ments while it is transported throughout the d isk (Dulle- 
mond & Dominik 2005 Birnstiel et al.||2010). A funda¬ 
mental limiting factor to the largest gram size is that of 
the radial drift problem. While solids move at Keple- 
rian velocities in the disk, the gaseous component moves 
at slightly sub-Keplerian speed due to gas pressure sup¬ 
port. This velocity difference induces a drag in the large 
particles, which end up losing angular momentum and 
drifting radially to smaller and smaller orbits, until they 
are accr eted onto the star and lost from the solid pop¬ 
ulation (Weidenschilling 1977). This barrier will limit 
the maximum gram size allowed in a disk, with mm and 
cm-sized particles at tens of AU from the star dri fting in¬ 
wards in timescale shorter than the disk lifetime (Youdin 


20101. Thus, measuring changes in the dust properties as 
a function of radius is essential to investigate the effects 
of this barrier in the grain size distribution of disks. 

Only recently have measurements reached the required 
sensitivity and a sufficient lever arm in wavelength 
to characterize radial variations of the dust properties 
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Guilloteau et al.||2011 Banzatti et al. 

20111 IFerez et al.| 

201211'Lrotta et al.||20l3| (Andrews et 


that spanned an order of magnitude, from sub-mm to 
cm wavelengths, observational constraints in dust opac¬ 
ity were obtained for the disk around the young star 
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j3 > 1.5 at more than ~80 AU, while a constant value of 
the dust opacity throughout the disk failed to reproduce 
these observations. Furthermore, these results seem to 
agree with a population of dust grains limited by radial 
drift. Here we explore further this issue, by studying 
how the dust properties vary with radius in the proto¬ 
planetary disks that surround the young stars CY Tau, 
located in the Taurus star fo rming region at a distance of 


140 pc (Torres et al. 2007), and DoAr 25, lo cated in the 
L1688 dark cloud at a distance of ~ 125 pc (Loinard et 


al. 2008 Mamajek 2008). DoAr 2f^ and CY Tau are 
pre-mam sequence stars of K5 anoMl spectral type, 


respectively. B oth stars are c 


uite young: CY Tau is 


2007 Guilloteau et al. 


2 — 3 Myr old (Bertout et al. 

20 14t, w hile DoAr 25 is about 
ar]|2009|, and both present a significant emission excess 


4 Myr old (Andrews et 


over the stellar phot osphere, from near-inf r ared to mil¬ 


li meter w avelengths (Olofsson et al. 


al. 2010). DoAr 25 has been imaged at su 


2009 


wa velengths w ith the Submillimeter Array (Andrews et 


al. 


McClure et 
D-millimeter 


2008 20091, and CY Tau was imaged at the 1.3 and 


Tg mm bands with the Plateau de Bure Interferometer 
( Guilloteau et al.||2011 ). 

in this paper, we triple the number of circumstellar 
disks for which our multiwavelength analysis has been 
used to measure any radial changes in the dust proper¬ 
ties, particularly in the dust opacity. The paper is struc¬ 
tured as follows: first, we describe the observations and 
data calibration procedures in Section and we present 
observational results from these data in Section 1 In 
Section |4] we discuss the model of the disk emission em¬ 
ployed in this paper, followed by Section where obser¬ 
vational constraints derived from this moaeling are pre¬ 
sented. Finally, in Section we discuss our results in the 
context of grain growth, comparing these observational 
constraints with the radial drift and fragmentation bar¬ 
riers for growth, as well as with other disks previously 
studied. Our findings and conclusions are presented in 
Section 0 


2. OBSERVATIONS 

We obtained multi-wavelength observations of the 
continuum emission from DoAr 25 and GY Tau with 
three interferometers: the Gombined Array for Re¬ 
search in Millimeter-wave Astronomy (CARMA), the 
Sub-Millimeter Array (SMA), and the Karl G. Jansky 
Very Large Array (VLA) as part of the Disks@EVLA 
collaboration. CY Tau was observed at four different 
wavelengths: 1.3 mm, 2.8 mm, 7.1 mm, and 5.0 cm, with 
the highest angular resolution being 0.07" (9.8 AU) at 
7.1 mm. DoAr 25 was observed at five different wave¬ 
lengths: 0.9 mm, 2.8 mm, 8.0 mm, 9.8 mm, and 5.0 cm, 
with the highest angular resolution being 0.10" (12.5 AU) 
at 8.0 mm. A summary of these observations and our 
data calibration procedure is presented below. 

2.1. CARMA Observations 

CY Tau and DoAr 25 observations at 2.8 mm 
(107 GHz), and GY Tau observations at 1.3 mm 
(230 GHz) were obtained with CARMA in the A, B, and 
C configurations, providing baseline lengths between 30- 
2000 m which cover spatial scales from 20" down to 0.3" 
at 2.8 mm, and between 9" down to 0.13" at 1.3 mm. 
A complete observing log of these observations can be 
found in Table [l] Double-sideband single-polarization 
receivers were tuned to a frequency of 107 GHz in the 
case of 2.8 mm observations, and 230 GHz in the case of 
1.3 mm observations. To optimize the continuum sensi¬ 
tivity, we configured all spectral windows in the correla¬ 
tor to the maximum possible bandwidth: 468.75-MHz- 
wide windows of 15 channels each, which provided 2.8 

also known as GY92 17, WSB 29, and WLY 1-34 
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TABLE 1 

Observing log for CARMA observations 


Target 

UT Date 

Configuration®' 

Phase Calibrator 


A = 1.3 mm 

CY Tau 

2007-NOV-07 

C, 0.35 km 

J0530-I-1331 


2009-Jan-04 

B, 1.0 km 

3C 111 


2009-Jan-05 

B, 1.0 km 

3C 111 


2010-Sep-16 

D, 0.25 km 

J0336-I-3218 


2010-Sep-21 

D, 0.25 km 

3C 111, J0336-I-3218 


2011-Dec-12 

A, 2.0 km 

J0336-I-3218 

A = 2.8 mm 

CY Tau 

2008-Dec-11 

B, 1.0 km 

3C 111 


2010-Feb-15 

A, 2.0 km 

3C 111,40336-1-3218 


2011-Dec-30 

B, 1.0 km 

3C 111 


2012-Feb-08 

C, 0.35 km 

3C 111 


2012-Feb-19 

C, 0.35 km 

3C 111 

DoAr 25 

2010-Jan-06 

B, 1.0 km 

J1625-2527 


2010-Jan-09 

B, 1.0 km 

J1625-2527 


2010-NOV-29 

A, 2.0 km 

J1625-2527 


2010-NOV-30 

A, 2.0 km 

J1625-2527 


2010-Dec-03 

A, 2.0 km 

J1625-2527 


2010-Dec-12 

A, 2.0 km 

J1625-2527 


2011-Apr-25 

C, 0.35 km 

J1625-2527 


2012-Jan-19 

C, 0.35 km 

J1625-2527 


^ Configuration name followed by the maximum baseline of each 
configuration. 


TABLE 2 

Observing log for VLA observations 


Target UT Date Configuration'^ Phase Calibrator 



A = 

7.1 mm (Q band) 


CY Tau 

2010- N0V-I3 

2011- Apr-05 

2012- Oct-22 
2012-Oct-26 
2012-Oct-27 
2012-Oct-28 

C, 3.4 km 

B, 11.1 km 

A, 36.4 km 

A, 36.4 km 

A, 36.4 km 

A, 36.4 km 

J0403-I-2600 

J0403-I-2600 

J0403-I-2600 

J0403-I-2600 

J0403-I-2600 

J0403-I-2600 

A = 8.0 and 9.8 mm (Ka band) 

DoAr 25 

2011-Jan-23 
2011-May-28 
2011-Jun-15 

CnB, 11.1 km 

BnA, 36.4 km 

A, 36.4 km 

J1625-2527 

J1625-2527 

J1625-2527 


A = 

5.0 cm (C band) 


CY Tau 
DoAr 25 

2011-Jul-23 

2011-Jul-14 

A, 36.4 km 

A, 36.4 km 

J0403-I-2600 

J1625-2527 


Configuration name followed by the maximum baseline of 
each configuration. 


GHz of bandwidth for observations before the correlator 
upgrade in December 2009, and 487.5-MHz-wide win¬ 
dows of 39 channels each, which provided 7.8 GHz of 
total bandwidth after the upgrade. 

The observing sequence interleaved observations of a 
complex gain calibrator (generally 3 minutes in length) 
with the science target. In C configuration the calibrator- 
target cycle time was 12-15 minutes, while in A and B 
configuration it was 5-10 minutes in order to track better 
the tropospheric phase fluctuations. This observing se¬ 
quence repeated throughout the track, usually from tar¬ 
get rise to set. In addition, a strong calibrator was ob¬ 


served to measure the complex bandpass, along with ei¬ 
ther a planet (Uranus, Mars, or Neptune) or a secondary 
flux density calibrator (either 3G273 or 3G84) that is 
monitored by the observatory. The estimated fractional 
uncertainty in the absolute flux density scale is ^ 15%. 

The GARMA data were calibrated using the Multi¬ 
channel Image Reconstruction, Ima ge Analysis and Dis - 


play (MIRIAD) software package (iSault et al.|fl995). 

Ma 


Each observation was calibrated separately. Maltunc- 
tioning antennas, receivers, and/or spectral windows 
were flagged, and updated antenna position corrections 
and line-length system corrections were applied. Obser¬ 
vations of DoAr 25 in the most extended (A and B) con¬ 
figurations used the GARMA Paired Antenna Calibra¬ 
tion System (C-PACS) to monitor the tropospheric delay 
fluctuations on eight of the outermost stations of the ar¬ 
ray, using adjacent 3.5 m antennas equipped with 1 cm 
receivers. Phase corrections derived from the C-PACS 
system are applied during post-proce ssing. A complete 
desc ription of C-PACS is p resented by Perez et al. (2010| 
and Zauderer et al. (2014). 


2.2. VLA Observations 
2.2.1. High Frequency 

CY Tau observations at 7.1 mm (42 GHz, Q-band) 
were made during A, B and C configurations of the VLA, 
while DoAr 25 observations at 8.0 and 9.8 mm (30.5 Ghz 
and 37.5 GHz, Ka-band) were made during A, BnA, and 
CnB configurationsj^providing baseline lengths for both 
targets between 60 m to 36 km. These baseline lengths 
correspond to spatial scale coverage from about 24" down 
to 0.04" at 7.1 mm, and from about 34" down to 0.06" 
at 9.8 mm. These observations, summarized in Table 
used dual-polarization receivers and two independently 
tunable basebands, where each baseband was configured 
to eight 128 MHz spectral windows of 64 channels, to pro¬ 
vide the maximum continuum bandwidth per baseband 
(1 GHz). For the Q-band observations, the two base¬ 
bands were centered at 7.1 mm (42.5 GHz) and 7.2 mm 
(41.5 GHz) providing 2 GHz total bandwidth at 7.14 mm, 
while for the Ka-band observations, the two 1 GHz base¬ 
bands were centered at 8.0 mm (37.5 GHz) and 9.8 mm 
(30.5 GHz). 


2.2.2. Low Frequency 

To distinguish the amount of free-free or non-thermal 
contamination present at millimeter and centimeter 
wavelengths, CYTau and DoAr 25 were observed at 5.0 
cm (6.0 GHz) in C-band. These low-frequency VLA ob¬ 
servations were obtained during A configuration, with 
dual-polarization receivers and two independently tun¬ 
able basebands, where each baseband was configured as 
for the high frequency observations with a total band¬ 
width per baseband of 1 GHz. The two basebands were 
centered at 6.2 cm (4.8 GHz) and 4.1 cm (7.3 GHz), in 
order to be able to infer the spectral slope of the emission 
at centimeter wavelengths. 

The hybrid BnA and CnB configurations comprise antennas 
on the east and west arms of the VLA moved to their locations for 
the more compact configuration, while the antennas on the north 
arm are moved to their locations for the more extended configu¬ 
ration, so that in projection, for southern sources, the synthesized 
beam is circularized. 
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Fig. 1. — Aperture synthesis images of the continuum emission towards the young star CYTau, observed at wavelengths of 1.3, 2.8, and 
7.1 mm. Each panel encompasses a 2.5” X 2.5” r^ion (350 AU in size at the adopted distance), with contours drawn at 3(t intervals, where 
cr is the RMS noise level in each map (see Table 


2.2.3. VLA data acquisition and calibration 

The observing sequence interleaved observations of a 
complex gain calibrator (45 sec to a few minutes in 
length) with science target observations, whose duration 
depended on the array configuration and observing fre¬ 
quency: target-calibrator cycle times in extended con¬ 
figurations were 1.5-5 minutes at high frequencies, 5-10 
minutes in compact configurations. The target-calibrator 
cycle time at C-band in A-configuration was 10-15 min¬ 
utes. A strong calibrator was used to determine the 
bandpass shape. The absolute flux density calibration 
was determined from observations of a primary flux den¬ 
sity calibrator (3C147 for CYTau, 3C286 for DoAr25). 
Models for the primary flux density calibrators are pro¬ 
vided by the observatory. The estimated uncertainty in 
the absolute flux density scale is 5% at C-band, and 10% 
at Ka- and Q-bands. 

The VLA data calibration was performed using the 
CASA software package and a modified version of the 
VLA calibration pipelinj^ As with the CARMA obser¬ 
vations, each observation was calibrated separately. Mal¬ 
functioning antennas, receivers, and/or spectral windows 
were flagged, as well as noisy channels at the edge of each 
spectral window. The first few seconds at the beginning 
of each target observation were also flagged, along with 
any radio frequency interference (especially an issue with 
the C-band data). Times of poor phase coherence, as 
measured on the gain calibrator scans, were also removed 
from the data. Because the data were acquired over 
an extended time period, corrections for source proper 
motion or other systematic position offsets were applied 
prior to combining. All VLA data presented here were 
obtained as part of the Disks@EVLA project (AC982). 


2.3. SMA Observations 

DoAr 25 observations at 0.9 mm (345 GHz) were ob¬ 
tained with the SMA between 2005 May and 2007 June; 
these observations have a lready been presented by |An- 
drews et al. (2008 2009). Three different array con- 
hgurations (C, E, and V) were used, providing baseline 
lengths between 8-500 m, corresponding to spatial scales 
from about 22" down to 0.3". Double-sideband receivers 
were tuned to a local oscillator frequency of 340.755 GHz. 
Each sideband was configured to have 24 partially over¬ 


processing/ pipeline 


https: / / science.nrao.edu / facilities/vla/data- 


lapping 104 MHz chunks, for a total continuum band¬ 
width of 4 GHz. The observing sequence interleaved 
observations of a complex gain calibrator with science 
target observations twice as long. The source-calibrator 
cycle was 8 minutes for the extended configuration, and 
^ 15-20 minutes for compact configurations. Bandpass 
and flux density calibrators were selected from different 
planets and satellites (Uranus, Jupiter, Saturn, Titan, 
Callisto), as well as strong quasars (3G454.3, 3C279), 
depending on their availability and array configuration. 
The estimated uncertainty in the absolute flux density 
scale is ^ 10%. The data were flagged and calibrated 
with the IDL-based MIR software package. 

2.4. Averaging of interferometric data 

For each telescope we verified that the calibration from 
multiple observations obtained in different array config¬ 
urations matched, to within the uncertainty in the ab¬ 
solute flux scale, by comparing the calibrated visibili¬ 
ties where they overlap in uv-sjoace. The observations 
at 0.9, 1.3, 2.8, 7.1, 8.0, and 9.8 mm were then aver¬ 
aged over the observing bandwidth, as follows. In the 
case of the 0.9, 1.3, 2.8 and 7 mm data, visibilities en¬ 
compass a narrow range of frequencies over the entire 
bandwidth: AA/A ^ 0.01 at 0.9 mm, AA/A ^ 0.04 at 
1.3 mm, AA/A ~ 0.08 at 2.8 mm, and AA/A ^ 0.05 at 
7.1 mm. These data were therefore averaged into a single 
wideband channel. In the case of Ka-band VLA obser¬ 
vations, the two 1-GHz wide basebands are separated by 
^ 7 GHz, so these data were averaged into two wide¬ 
band channels: one centered at 30.5 GHz (A = 9.8 mm), 
the other centered at 37.5 GHz (A = 8.0 mm). Gali- 
brated observations from different array configurations 
were combined to generate a single visibility file for each 
wavelength analyzed, in order to compare these interfer¬ 
ometric observations with physical disk models. 

3. OBSERVATIONAL RESULTS 
3.1. CY Tau 

Synthesized maps at 1.3, 2.8, and 7.1 mm of the con¬ 
tinuum emission from CYTau are presented in Figure 
Each map extends 2.5" x 2.5", corresponding to 350 AU 
at the adopted distance. At wavelengths of 1.3 mm, 
7.1 mm, and 5.0 cm, natural weighting was used to max¬ 
imize the sensitivity. At a wavelength of 2.8 mm, a 
Briggs weighting scheme with a robust parameter of 0.5 
was used in CASA, to optimize a combination of resolu- 
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Fig. 2.— Aperture synthesis images of the continuum emission towards the young star DoAr 25, observed at wavelengths of 0.9, 2.8, 8.0, 
and 9.8 mm. Each panel encompasses a 2.8^' x 2 . 1 " region (350 AU in size at the adopted distance), with contours drawn at 3cr intervals, 
where cr is the RMS noise level in each map (see Table]^. 


tion and sensitivity. The resulting image properties and 
source photometry can be found in Table from these 
measurements we infer an spectral index for^Y Tau from 
1.3 to 7.1 mm of a = 2.6. 

The observations at 5 cm were used to estimate the 
contribution from processes other than thermal dust 
emission at 7.1 mm (e.g., chromospheric activity or ther¬ 
mal bremsstrahlung from phot oevaporative disk winds 
driven by the cen tral protostar; |Mundy et al.|1993{|Pas-| 
cucci et d^|2011 ), which needs to be taken into account 


when modeling the dust emission from the disk. How¬ 
ever, GY Tau is not detected at 5 cm, thus we take the 3 (t 
upper limit of 20 /rJy and estimate the maximum possi¬ 
ble contamination at 7.1 mm by assuming optically-thick 
free-free emission with spherical symmetry, for which the 
emission as a function of fr equency, is proportional to 
Si, oc (Reynolds 1986). The maximum contamina¬ 
tion at 7.1 mm thus corresponds to ~ 60 ^Jy, which 
represents only 4% of the total emission at this wave¬ 
length. We note that the probable source of any such 
contamination will arise from very near the protostar, 
and would appear as an unresolved point source in the 
7.1 mm data. Based on the 7.1 mm data themselves, us¬ 
ing Mu-distances > 10 km (corresponding to spatial scales 
smaller than ^ 0.14"), we constrain any unresolved point 
source to have a flux density of 74±9 /iJy (some of which 
could be dust emission), consistent with our estimate of 
the potential contamination extrapolated from 5 cm. 


3.2. DoAr 25 

Synthesized maps at 0.9, 2.8, 8.0 and 9.8 mm of the 
continuum emission from CY Tau are presented in Fig¬ 
ure Each map extends 2.8" x 2.8", corresponding to 
350WU at the adopted distance. Natural weighting was 
used in the imaging at wavelengths of 8.0 mm, 9.8 mm, 
and 5.0 cm, to maximize the sensitivity. Briggs weight¬ 
ing with robust parameters of 0.7 and 0.3 were used for 
imaging at wavelengths of 0.9 and 2.8 mm, respectively, 
to optimize the resolution and sensitivity at these wave¬ 
lengths. The resulting image properties and source pho¬ 
tometry can be found in Table from these measure¬ 
ments we infer an spectral index for DoAr 25 from 0.9 to 
9.8 mm of a = 2.8. 

As for CY Tau, the observations at 5 cm from DoAr 25 
were used to estimate the contribution from other than 
thermal dust emission at 8.0 and 9.8 mm. In the case 
of DoAr 25, A = 5 cm emission is detected at a level of 


6(7 coincident with the star, with an integrated flux den¬ 
sity of 50 ± 13 ^Jy. Assuming (x ® results in an 
estimated contamination of 150 ±40 /iJy at 8.0 mm, and 
133 ± 35 /LtJy at 9.8 mm, which corresponds to 17% and 
24% of the integrated emission at these wavelengths. A 
lower level of contamination is derived based only in the 
8.0 and 9.8 mm data with uiedistances > 10 km (cor¬ 
responding to spatial scales smaller than ~ 0.17"): an 
unresolved point source flux density of 68 ± 24 /xJy is 
present in the 8.0 mm observations, while a flux den¬ 
sity of 65 ± 16 /rJy is constrained at 9.8 mm. These 
estimates are consistent (within the errorbars) with the 
extrapolation from 5 cm. However, we adopt the most 
conservative value (the highest possible contamination) 
in our analysis. 


3.3. Brightness temperature and optical depth of the 
emission 

The brightness temperature of the emission at mm and 
cm wavelengths can be useful to discriminate between op¬ 
tically thin and optically thick emission. For a medium 
at a physical temperature T with no background radia¬ 
tion and an optical depth r, the brightness temperature 
(Tb) of the emission at a particular wavelength A is re¬ 
lated to r as: Tb{\) = T{1 — (see, e.g., |Wilson 

et al. 2009). Thus, in the optically thick limit (r ^ 1) 
the brightiiess temperature directly traces the tempera¬ 
ture of the medium and is independent of wavelength: 
T'b(A) = T, while in the optically thin limit (r ^ 1) 
the brightness temperature will be much lower than the 
physical temperature of the medium: Tb(A) = t\T and 
will vary with wavelength according to the wavelength 
dependence of t\. 

The data presented in §3.1 and §3.2 include high fre¬ 
quency observations and emission from the cold outer 
disk of CY Tau and DoAr 25, thus the Rayleigh-Jeans 
limit of hv kT might not be appropriate even at cm 
wavelengths. For this reason we calculate the brightness 
temperature of the emission using the full Planck func¬ 
tion without this approximation. Figuresj^and ^present 
azimuthally-averaged radial Tb profiles for botn disks at 
each observed wavelength, with the shaded region rep¬ 
resenting the ItT constraint derived from the statistical 
uncertainty in the images. For a proper comparison, we 
convolved all observations to the same angular resolution 
(i.e., the lowest spatial resolution, which corresponds to 
the resolution of the 2.8 mm data). In addition, we sub- 
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TABLE 3 

CYTau imaging and photometry results 





Image properties 


Source photometry^ 

Telescope 

A 

[mm] 

Image rms noise 
[/ijy/beam] 

Synthesized beam 

Beam P.A. 

n 

Flux density 
[mjy] 

CARMA 

1.3 

580 

0.29" X 0.24" 

83.3 

119 ± 18 

CARMA 

2.8 

280 

0.40" X 0.35" 

-87.1 

26.1 ± 3.9 

VLA 

7.1 

7.8 

0.07" X 0.07" 

-49.5 

1.76 ±0.27 

VLA 

50 

7.0 

0.62" X 0.34" 

-63.7 

< 0.020 (3cr) 


^ Source photometry obtained from the measured flux density of CY Tau limited to un-distances 
between 0-40 kA, corresponding to the shortest un-spacings we have sampled, except for A = 50 mm 
which is not detected and thus we can only provide an upper limit. Note that the errors quoted 
include the uncertainty in the absolute flux density scale. 


TABLE 4 

DoAr25 imaging and photometry results 


Telescope 

A 

[mm] 


Image properties 


Source photometry^ 

Image rms noise 
[/ijy/beam] 

Synthesized beam 

Beam P.A. 

[°] 

Flux density 
[mJy] 

SMA 

0.9 

3500 

0.48" X 0.35" 

14.4 

515 ± 52 

CARMA 

2.8 

270 

0.64" X 0.33" 

0.9 

28.1 ±4.2 

VLA 

8.0 

11.0 

0.15" X 0.10" 

-13.2 

1.17 ± 0.13 

VLA 

9.8 

8.0 

0.18" X 0.12" 

-13.3 

0.66 ±0.07 

VLA 

50 

5.9 

0.68" X 0.34" 

-8.4 

0.051 ± 0.012 


^ Source photometry obtained from the measured flux density of CY Tau limited to un-distances 
between 0-40 kA, corresponding to the shortest un-spacings we have sampled except for A = 50 mm, 
where the source is weak and unresolved and the fit is performed in the image domain. Note that 
the errors quoted include the uncertainty in the absolute flux density scale. 


tracted an unresolved point-source component for the 
long-wavelength observations at 7.1 mm (CYTau) and 
at 8.0 and 9.8 mm (DoAr25), whose flux density is given 
by the amount of contamination from non-dust emission 
described in §3.1 and §3.2. 

The observed brightness temperature profiles in Fig¬ 
ures and 1^ decline with increasing distance from the 
star in both disks, becoming successively fainter as the 
wavelength of the emission is increased. If the observed 
emission were optically thick, the Tg profiles would di¬ 
rectly trace the physical temperature of the disk. For 
a vertically isothermal disk in this limit, the observed 
brightness temperature should then be the same at all 
wavelengths. For both DoAr25 and CYTau we observe 
a successively fainter Tb profile with increasing wave¬ 
length, suggesting that for all but the shortest wave¬ 
length, the observed emission is consistent with being 
optically thin. Furthermore, in §5 we find that, for a 
flared disk in hydrostatic equilibrium, the midplane tem¬ 
perature is a factor of several higher than the observed 
brightness temperature of CY Tau and DoAr 25 at all 
disk radii and at all wavelengths, thus indicating that 
the observed emission is optically thin rather than opti¬ 
cally thick (see section §4 for details of the modeling and 
§5.1 for the resulting midplane temperature profiles). 

A colder midplane temperature could be found in the 
extreme case of a thin flat disk in the presence of no 
interstellar radiation. This configuration results in the 
coldest midplane disk possible, since the opening angle of 
a thin disk is much smaller than for a flared disk, making 
the stellar heating less efficient. Note that a more settled 
disk interpretation is favored in the analysis of unresolved 
observations of CY Tau, as the SED from mid-infrared 


to far-infr ared wavelengths is q uite steep for this object 
(see, e.g., |Robitaille et al.|2007|. In the case of DoAr 25, 
Spitzer IRS spectra reveal a flat spectrum betw een 10-30 
microns consistent with a flared disk structure (|01ofsson|| 
et al.|2009( |McClure et al.|2010|). Thus, we only consider 
the thin disk case for CY "i'au. We find a midplane tem- 


midplane 

perature of ~ 14 K at 10 AU, decreasing to a lower limit 
defined by interstellar radiation and cosmic ray heating 
(assumed to be ^ 8 K) at radii of 30 AU and greater. 
Only for this extreme case is the observed Tb profile of 
CY Tau consistent with the physical temperature of the 
disk midplane, and this is only for the shortest wave¬ 
length observations at 1.3 mm at radii ^30 to 50 AU. 
At longer wavelengths (and at larger radii at 1.3 mm) 
the brightness temperature is still below this cold mid¬ 
plane case, which indicates that the optical depth of the 
emission is less than 1, even in this extreme case. Con¬ 
sequently, we will assume in the following analysis that 
the emission at 1.3 mm from CYTau is optically thin, 
although there is a possibility of higher optical depth at 
1.3 mm if the disk structure corresponds to a flat thin 
disk for this object. 


3.4. De-projected visibility profiles 

In Figures and we present the real and imaginary 
part of the visibility as a function of wn-distance (here¬ 
after referred to as visibility profiles), for both CYTau 
and DoAr 25, at each of the observed wavelengths. Each 
visibility was deprojected by the known inclination and 
position angle of the disk (see Section]^ for details of the 
known disk geometry), before averaging into uiubins with 
a width of 40kA. To compare the different wavelengths, 
each visibility profile is normalized by the measured flux 
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CY Tau 



Fig. 3. — Observed brightness temperature of the emission from 
our multi-wavelength observations of CYTau, as a function of the 
radial distance to the central star. These profiles are being com¬ 
pared at the same angular resolution, corresponding to the 2.8 mm 
observations (0.4^^, 56 AU). 


DoAr 25 



Fig. 4.— Observed brightness temperature of the emission from 
our multi-wavelength observations of DoAr 25, as a function of 
the radial distance to the central star. These profiles are being 
compared at the same angular resolution, corresponding to the 
2.8 mm observations (0.64^^, 80 AU). 


density in the first uv-hin^ between 0-40 kA. In Fig- 
ures and the black dotted line at a constant value of 
Re = 1.0 represents the observed visibility profile for an 
unresolved point source. Once a disk is resolved by inter¬ 
ferometric observations, its visibility profile will decline 
from this reference line. Thus, the significant decline of 
the real part of the visibility profile at each wavelength 
demonstrates that the emission from both DoAr 25 and 
CYTau is resolved, at all the observed wavelengths. Fur¬ 
thermore, a compact source will have a shallower decline 
with M?;-distance than an extended source, whose decline 
will be quite steep. Thus, since the decline of the short- 
wavelength visibility profile (0.9, 1.3, and/or 2.8 mm) is 
steeper than the long-wavelength visibility profile (7.1, 
8.0 and/or 9.8 mm), the disk emission observed at short- 
wavelengths is more extended than the disk emission ob¬ 
served at long-wavelengths, which has to be more com¬ 


pact. Hence, the differing visibility profiles demonstrate 
a wavelength-dependent structure in both the CY Tau 
and DoAr 25 disks, s imilar to that fou nd for the disks 
surrounding AS 209 (Perez et al. 2012), CQ Tau (Ban- 
zatti et al. 2011), and several other stars in the Taurus- 


Auriga star-forming region (Guilloteau et al. 2011). In 
the next section, we explaiii the wavelength-dependent 
structure of the emission from CY Tau and DoAr 25 as 
radial variations of the dust properties across the circum- 
stellar disk. 


4. MODELING OF THE DISK EMISSION 

Observations of CY Tau and DoAr 25 were analyzed 
using a disk emission model that reproduces the ra¬ 
dial brightness distribution at mm and cm wavelengths, 
equivalent to the disk model employed in the analysis of 
similar obser vations of the disk surrounding the young 
star AS 209 (iPerez et al.||2012|). At A = 7.1, 8.0, and 
9.8 mm an additional point source at the center of the 
disk is also included, to represent the potential contami¬ 
nation from sources othe r th an th ermal dust emission, as 
described in Sections |01 and [3)^ We adopt the two-layer 
dis k approximation for the disk structure, first presented 
by Chiang & Goldreich (1997). In this model, a flared, 
passively-heated disk is irradiated by the central star, 
with its structure characterized by two components: a 
disk surface layer that absorbs the stellar radiation, and 
a disk interior, which is opaque to the stellar photons. 
Dust in the surface layer will absorb the short-wavelength 
stellar emission and re-radiate it at longer wavelengths. 
As the disk interior absorbs the reprocessed emission, the 
inner regions of the disk can be efficiently heated. We 
use the implementa tion of this r adiati ve transfer prob¬ 
lem as described by Isella et al. (2009), where the tem¬ 
perature structure of the two-layer m odel is computed 
with t he iterative method presented by |Dullemond et al.| 
(2001). The assumptions implicit in the two-layer disk 
model break down for the far outer regions of the disk, 
where the surface density of material drops significantly 
and even the disk interior becomes optically thin to the 
stellar radiation. When this condition is met, we assume 
the same temperature for both the surface layer and disk 
interior. The disk temperature smoothly decreases as a 
power-law function of radius down to a minimum of 10 K 
in the outer disk, where the dust is heated by the inter¬ 
stellar radiation field. 

Por the mass surface density structure, E(i?), we adopt 
the similarity solution for a viscously-accreting Keplerian 
disk subject to the gravity of a massive central object 
(Lynden-Bell & Pringle 1974), which can be described 
by a power law combined with an exponential taper at 
large radii. If the disk viscosity, Vy, is assumed to be a 
power-law with radius {vy oc R'^), then the s olution for 
the disk su rface density is time-independent (Hartmann 
et al.|1998 ). We adopt the sam e parameterization ot th e 
similarity solution presented by Guilloteau et al. (2011): 


S(i?) = Eo ( ^ 
, rto 


X exp — 


R 


2 - 7 ^ 


( 1 ) 


where Eq corresponds to the surface density at radius 
i?o, 7 is the power-law exponent on the radial depen¬ 
dence of the viscosity, and Rc is the characteristic seal- 

















































Perez et al. 



0 200 400 600 800 1000 1200 1400 
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Fig. 5. — Normalized real and imaginary part of the correlated 
emission from CY Tan as a function of spatial frequency {uv- 
distance). The visibilities have been de-projected using a position 
angle of 63 ° and inclination of 28° , derived from molecular line ob¬ 
servations jGuilloteau et al.|2011|. Each bin has a width of 40 kA, 
and has been normalized by tne measured flux density at 0-40 kA. 
Filled circles and error bars of different color correspond to corre¬ 
lated real and imaginary part of the emission observed at different 
wavelengths (1.3 mm: blue circle, 2.8 mm: gray triangle, 7.1 mm: 
black square). The continuous lines correspond to the best-fit disk 
emission model at each wavelength (see modeling in §5.1). 



uv~distance [kA] 

Fig. 6 . — Normalized real and imaginary part of the corre¬ 
lated emission from DoAr 25 as a function of spatial frequency 
(ui;-distance). The visibilities have been de-projected assuming 
a position angle of 109° and inclination of 62°, as derived from 
the modeling of the shortest-wavelength emission. Each bin has a 
width of 40 kA, and has been normalized by the measured flux den¬ 
sity at 0-40 kA. Filled circles and error bars of different color cor¬ 
respond to the correlated real and imaginary parts of the emission 
observed at different wavelengths (0.9 mm: blue circle, 2.8 mm: 
gray triangle, 8.0 mm: black square, 9.8 mm: cyan diamond). The 
continuous lines correspond to the best-fit disk emission model at 
each wavelength (see modeling in §5.1). 


ing radiuEl This prescription for S(-R) behaves as a 
power law for small radii, while at large radii it decreases 
smoothly in an exponential fashion. 

The inner radius of the disk, Rin, is dehned as the 
radius at which the dust sublimates: 


= 0.07 AU ^ 


(a 


1/2 / rji 

I ^ subl 


Vl500 K 


-2 


( 2 ) 


which corresponds to ^ 0.04 AU for CY Tau and ^ 
0.06 AU for DoAr 25, for a dust sublimation tempera¬ 
ture of Tsubi = 1500 K. Given that the angular resolution 
of our observations is more than an order of magnitude 
greater than Rim the value of this parameter does not 
affect our constraints on the disk structure. 


As in the analysis presented by Perez et al. (2012), 
we compute the magnitude of the dust opacity, us¬ 
ing Mie theory. We adopt a grain population of com¬ 
pact spherical grains larger than Umm = 0.01 /xm, in 
a power-law distribution of grain sizes: n(a) oc a~^ for 


The scaling radius Rq is directly related to the location in 
the disk where the mass accretion rate changes sign, the transition 
/ \i/( 2 — r) 

radius: Rt = Rc ( 2 ( 2 - 7 ) / ' transition radius defines a 

boundary for accretion or expansion, since for R < Rt mass flows 
inward while for R > Rt mass flows outward. Note that in the 
context of viscous disk evolution, the value of the transition radius 
will increase with time due to conservation of angular momentum: 
as matter is accreted onto the star the disk must expand to con¬ 
serve total angular momentum. We include its definition here since 
different authors will either use Rq or Rt in their preferred pre¬ 
scription for the disk surface density. 


o-min < a < Umax- Grains in the disk interior are as¬ 
sumed to be composed of silicates, organic materials, and 
water ice, while grains in the disk surface are assumed to 
be depleted of ice. We obtain optical cons tants for these 
grain materia l s from Semenov et al. (2003) (for silic ates) 


Zubko et al. 


ren 


( 1996[) (tor amorphous carbon), and War- 


]iyS4r(±or water ice). Rec ent studies have rev ised 
the solar abundance of oxygen ( Asplund et al.|[200^ , so 
although we assume the fractional abund ances ot these 
materials as given by Pollack et al. (1994), we reallocate 
the fractions of these astrophysical grams to account for 
the increased oxygen abundance. This results in grains 
in the disk interior which are composed of 12% silicates, 
44% organics and 44% water ice, while in the disk sur¬ 
face (where water ice may not be present) the fractional 
abundances correspond to 21% silicates and 79% organ¬ 
ics. Except for the normalization of the surface density 
at radius Rg (Eq); the effect of dust grain composition in 
the disk modeling will produce a minimal effect in the de¬ 
rived parameters of the disk structure (well with in their 
uncertainties, as shown by Isella et al. 2010). More 
specifically, although the dust composition does affect 
the absolute value of opacity, the fact that the dust emis¬ 
sion is optically thin at long wavelengths (§5.3) means 
that the derivation of the dust opacity spectral index /3 
is insensitive to the assumed absolute opacity value and 
hence dust composition. This follows since the only fre¬ 
quency dependence when the emission is optically thin 
corresponds to oc * B^{T). 
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Finally, the last parameters that determine the ob¬ 
served disk emission describe the geometry of the disk 
in the plane of the sky: inclination (i) and position an¬ 
gle (PA, measured from North to East). In the case 
of GY Tau, the disk geometry has been characterized 
from molecular line observation s of CO at 0.5" resolu - 
tion: i = 28°±5°, PA= 63° ±1° ( [Guilloteau et al.|201lD , 
which is what we adopted for this study (note that these 
constraints are consistent with those derived from CN 
observations at lower spatial resolution: i = 24.0° ±2.4°, 
PA= 62.5° ± 1.8°;|Guilloteau et al.||2014|). For DoAr 25, 


there are no molecular line observations in the literature. 
We therefore constrained its disk geometry from our 
modeling of the dust continuum observations at 0.9 mm, 
where the signal-to-noise ratio is the highest, finding 
i = 62°±1° and PA= 109° ±2° (c onsistent with previous 
modeling by | Andrews et al.||2008l |2009|). 

We find the best-ht model to a single wavelength obser¬ 
vation through minimization usin g a Markov-Chain 
Monte Carlo (MCMC) procedure (see Isella et al. ]200^ . 
The probability distribution is sampled by varying 
the free parameters that define the surface density pro¬ 
file (Sq, Rc,x), since the inclination and position an¬ 
gle of the disk are fixed to the values reported above, 
and where Rq has been fixed to a value of 20 AU. A set 
of these parameters defines a model of the disk bright¬ 
ness distribution, with the addition of a point source at 
the disk center whose flux density corresponds to the es¬ 
timated level of contamination from non-dust thermal 
emission, as inferred in Section!^ From this model emis¬ 
sion, we produce an image of the disk making sure that 
the significant spatial scales are covered. We take the 
Fourier transform of this image and obtain model visi¬ 
bilities sampled at the same locations in the Mw-domain 
as the original observations. We compare model and ob¬ 
served visibilities by means of the and construct 16 
separate MCMC chains that sample the parameter space 
using a Metropolis-Hastings algorithm with a Gibbs sam¬ 
pler. 

These chains all converged to the best-fit model, which 
corresponds to the one that minimizes the x^ and hence 
is the model that best reproduces our observations. We 
note that this MCMC algorithm samples the posterior 
probability distribution function (PDF) of the param¬ 
eters when chains have converged and reached equilib¬ 
rium. Thus, to find best-fit values and confidence inter¬ 
vals for the parameters in our model we adopt a Bayesian 
approach, marginalizing the resulting PDF over all but 
one parameter to obtain the probability distribution of 
the parameter of interest. From each parameter PDF we 
find the region that contains 68.3%(ltT), 95.5%(2 (t), and 
99.7%(3 (t) of all samples at equal probability to constrain 
each of the aforementioned parameters. 

5. MODELING RESULTS 


To model our CY Tau and DoAr 25 observations, we 
first select the maximum grain size (umax) and grain size 
distribution slope (g) that best reproduces the unresolved 
SED from sub-mm to cm wavelengths. For CY Tau, these 
correspond to q = 3.5 and a^ax = 2.5 mm, which results 
in an opacity spectral slope P = 0.79 between 1.3 and 
7.1 mm. For DoAr 25, the long-wavelength SED can be 
reproduced with q = 3.5 and Umax = 1-5 mm, resulting 
in an opacity spectral slope P = 0.91 between 0.9 and 


TABLE 5 

Best-fit model parameters and lu 

CONSTRAINTS FOR CY TAU 


A 

Rc 

7 

So 

X?ed 

[mm] 

[AU] 


[gm cm“^] 


1.3 

62.5ll;« 

U.UD_o 06 

^•oO_o,o7 

1.15 

2.8 

64.411? 

U.lD_o 07 

Q OQ + O.ll 
'^•'^^-0.12 

1.14 

7.1 

36.0l?;g 

f) kc:+0-06 

U.OO_o 08 

6 1 

1.14 


TABLE 6 

Best-fit model parameters and la constraints 
EOR DoAr 25 


A 

Rc 

7 

So 

X?ed 

[mm] 

[AU] 


[gm cm~^] 


0.9 

123.21?® 

—0 

^•-*-^-0.05 

^•o^_0.03 

1.08 

2.8 

103.5I4? 

0-34l?-°? 

01 +0.05 

-L-^-L_o.06 

1.07 

8.0 

IO6.0I??? 

0.83l?;f3 

1 QO + O-ll 
-*^•^^-0.08 

1.11 

9.8 

45.51?;? 

O.32I?;?? 

^•''^-0.28 

1.12 


8.0 mm. Note that while these values of Umax and /3 are 
consistent with the SEDs, they are relatively poorly con¬ 
strained by the curvature of the long wavelength emis¬ 
sion. 

In order to minimize the number of free parameters 
in this model, we assume a radially constant dust opac¬ 
ity throughout the disk and fit each wavelength sepa¬ 
rately. Our model then constrains the disk surface den¬ 
sity (S(i?)) and temperature profile {T{R)) at each wave¬ 
length; these constraints will appear to be different at dif¬ 
ferent wavelengths if the assumption of a constant dust 
opacity is not satisfied. The same approach has been 
successfully employed to measure radial variations of the 
dust opacity in the disks surrounding RY Tau, D C Tau, 
and AS 209 rtlsella et al.||2010t |Perez et al.||2012|); in the 
following sections we present our results tor the circum- 
stellar disks of CY Tau and DoAr 25. 


5.1. CY Tau and DoAr 25 disk structure constraints 

The best-fit and Icr constraints for each model parame¬ 
ter (Yq, Rc,x), obtained from their marginalized proba¬ 
bility distribution function, are presented in Tables and 
m To demonstrate that these best-fit models are sensi¬ 
ble representations of the observations, we imaged the 
model and residual visibilities (calculated by subtract¬ 
ing the model visibilites from the data) in the same way 
as the original observations. The resulting images are 
shown in Figure and Figure for CYTau and DoAr 25 
respectively. Additionally, for each of the best-fit models, 
we computed de-projected visibility profiles (solid lines 
in Figures and 1^. The adopted model is a good rep¬ 
resentation of the data at each separate wavelength, as 
no significant residual emission can be seen in the resid¬ 
ual maps and the visibility profile for our best-fit model 
agrees with the observations at all wavelengths. 

From each separate model fitting we characterize the 
disk surface density and temperature profile (E(R) and 
T(i?)) at each wavelength. Figures [9 and 10 present the 
best-fit and 3cr constraints on the disk temperature, sur¬ 
face density, and optical depth = k\Y{R)), obtained 
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Fig. 7. — Modeling of the 1.3 mm (top), 2.8 mm (middle), and 
7.1 mm (bottom) continuum emission towards CYTau. For each 
row, the data are shown on the left, the best-fit model disk is shown 
on the middle panel, and the residual (obtained by subtracting the 
best-fit model from the data in the Fourier domain) is on the right. 
The imaging parameters are as described in Section |3.1| Contours 
start at —3^ stepping by 3a, where a is the RMS noise in each 
map (Table 1^. 
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Fig. 8 .— Modeling of the 0.9 mm (top), 2.8 mm (2nd from top), 
8.0 and 9.8 mm (bottom two rows) continuum emission towards 
DoAr 25. For each row, the data are shown on the left, the best-fit 
model disk is shown on the middle panel, and the residual (obtained 
by subtracting the best-fit model from the data in the Fourier do¬ 
main) i s on the right. The imaging parameters are as described in 
Section |3.2| Contours start at —3f7, stepping by 3 (t, where a is the 
RMS noise in each map (Table [^. 


from modeling each wavelength independently. The tem¬ 
perature prohles inferred from observations at different 
wavelengths are very similar (this is expected, since dust 
of different sizes in the disk should have reached radiative 
equilibrium and be at the same temperature). In partic¬ 
ular, the best-fit temperature at each wavelength differs 
from the mean temperature profile by at most ^ 3 K for 
CY Tau and ~ 5 K for DoAr 25. However, the surface 
density profiles that were inferred separately from each 
observation (middle panels of Figures |9 and[l^, are dif¬ 
ferent at each observed wavelength. Tnis is clearly not 
physical (the dust emission arises from the same disk at 
all wavelengths!), and cannot be reconciled through a 
global change in the dust opacity. We conclude that we 
need to consider a change in the dust opacity as a func¬ 
tion of location (specifically, radius) in the disk, in order 
to explain the observed dust emission. This is explored 
further in §5.2| below. Finally, the right panels of Figures 

t and To] show the inferred optical depth, demonstrating 
at thedust emission is optically thin at all wavelengths 
for i? > 15 AU (well within our spatial resolution at all 
wavelengths for both disks). Furthermore, the emission 
would be optically thin even if the mid-plane dust tem¬ 
perature were lower by a factor of a few. 


5.2. Radial variations oj the dust opacity 

Since the dust emission is in the optically thin regime 
for both CYTau and DoAr25 (right panel of Figures]^ 
and 10), the observed emission will directly depend on 


three physical quantities: the dust opacity, the dust 
mass surface density, and the temperature, such that 
Si, cx k\T,B,j{T). Thus, our observations provide a con¬ 
straint for the product Ky x E x B,y{T), where the opac¬ 
ity, surface density, and temperature may vary with ra¬ 
dius, but only the dust opacity is a function of the ob¬ 
served wavelength. Our model fitting presented above - 
which assumes a constant dust opacity with radius and 
where each wavelength is fitted separately - results in a 
wavelength-dependent E(i?), signifying that the assump¬ 
tion of a radially constant dust opacity is not warranted 
by the data. To reconcile these wavelength-dependent 
T,\{R) and Ty(i?) (Figures and [^, as well as the 
wavelength-dependent disk structure seen in the nor¬ 
malized visibility profiles (Figures and we require 
a change in the dust properties as a function of radial 
distance from the central star, for both the CY Tau and 
DoAr 25 disks. 

We focus on the interpretation that the spectral index 
of the dust opacity, (3, is changing with radius, which 



unknown) disk structure: 


nx Ey(i?) B,{Tx{R)) = kx{R) E(i?) B,{T{R)) (3) 


where kx{R), T,{R), and T{R) are the true disk phys¬ 
ical quantities, and the left-hand side of this equation 
corresponds to our multi-wavelength constraints for the 
disk structure. Thus, two different wavelengths, Ai and 
A 2 (and corresponding frequencies, ui and 1 / 2 ), can be 
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CYTau 


CY Tau 


CY Tau 





Fig. 9.— Left and middle: Tx{R) and inferred from separate modeling of our multi-wavelength observations of CYTau, assuming a 

constant dust opacity with radius. Right: Optical depth tx(R) = i^xX T,x{R) inferred from separate modeling of CYTau multi-wavelength 
observations, assuming a radially constant kx- Colored regions: 3cr confidence interval constrained by our observations; continuous line: 
best-fit model; dashed line on left panel: mean temperature profile. The different S(H) and T(R) profiles for each wavelength imply a 
varying dust opacity with radius, and because of this, none of them is the true surface density and temperature profile of the disk. 



Fig. 10.— Left and middle: Tx{R) and T,x{R) inferred from separate modeling of our multi-wavelength observations of DoAr 25, assuming 
a constant dust opacity with radius. Right: Optical depth tx {R) = k-xX Tix{R) inferred from separate modeling of DoAr 25 multi-wavelength 
observations, assuming a radially constant kx- Colored regions: 3cr confidence interval constrained by our observation; continuous line: 
best-fit model; dashed line on left panel: mean temperature profile. The different S(H) and T(R) profiles for each wavelength imply a 
varying dust opacity with radius, and because of this, none of them is the true surface density and temperature profile of the disk. 


related by 


-/3(fl) 


-Pc Sai(-R) 


(R)) 


/ Bu^ {Tx 2 i^)) 


( 4 ) 


assuming that at long wavelengths the dust opacity fol¬ 
lows a power-law behavior {k\ oc A“^), and Pc corre¬ 
sponds to the spectral slope of the assumed constant 
dust opacity with radius, derived from the integrated 
SED used in the initial models. A useful prescription 
to infer radial variations of the opacity spectral slope 
(A/3(i?) = P{R) — Pc) in dual-wavelength observations 
is de rived from this las t equation (see also equation (9) 
from Isella et al. (2010)): 


A/3(i?) = - log 


-1 


Al 

A2 


X log 




( 5 ) 

where the dust surface density and temperature depend 
on the distance from the star R. Equation [^implies that 
in logarithmic space AP{R) is the slope of a straight line 
at points {a; = log(A),j/ = \og{i:xB^{Tx))/B^{T)}. 

Note that if the emission is in the Rayleigh-Jeans (R-J) 
domain: B^{T) oc T, and thus the true dust temperature 


(T{R)) will not be needed to derive AP{R). However, 
at cold temperatures in the outer disk and/or at short 
wavelengths, the R-J assumption breaks down and w e re- 
quire an estimate of T{R). As mentioned in Section 5.1 


the temperature constraints from these multi-wavelength 
data are quite similar, thus we assume the true dust tem¬ 
perature to be the mean TAR), sho wn as a dashed line 
on the left panels of Eigures[^and Note that a change 
in the estimate of T{R) will have a minimal impact in the 
AP{R) constraint, as T{R) only appears in the ratio of 
Planck functions evaluated at two different wavelengths 
in equation (and close to the R-J limit the dependence 


with T{R) is negligible). Nevertheless, we tested two dif¬ 
ferent assumptions where T{R) was equal to twice and 
half the mean of Tx{R). These substantial changes of 
x2 only changed our constraints in AP{R) by ~ 10% or 
less throughout the disk. 

Since the result of the MCMC algorithm is a sampled 
posterior PDF for the parameter space, we can construct 
a PDF for the product Ea x B^{Tx)/B^{T) at each radius 
and for each wavelength that has been modeled. To de¬ 
termine AP{R) we use a Monte Carlo simulation, where 
the PDF of Ea x B^{Tx)/B,^{T) is randomly sampled 
at each wavelength; a random sampling of AP{R) will 
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CY Tau 



Fig. 11. — Spectral index of the dust opacity (^) as a function of 
radius, inferred from multi-wavelength observations of the CY Tau 
disk. Black line: best-fit /3(il); colored areas: confidence interval 
constrained by our observations. The vertical dashed-lines repre¬ 
sent the smallest spatial scale probed at each wavelength, the arrow 
indicates that at 7.1 mm the smallest spatial scale is outside of the 
range plotted in this figure, corresponding to 42 mas or 6 AU. The 
errorbar at the bottom right indicates the systematic uncertainty 
on resulting from absolute flux scale uncertainty. 


DoAr 25 



Fig. 12.— Spectral index of the dust opacity (/3) as a function of 
radius, inferred from multi-wavelength observations of the DoAr 25 
disk. Black line: best-fit /3(/?); colored areas: confidence interval 
constrained by our observations. The vertical dashed-lines repre¬ 
sent the smallest spatial scale probed at each wavelength, the arrow 
indicates that at 8.0 and 9.8 mm the smallest spatial scale is out¬ 
side of the range plotted in this figure, corresponding to 48 mas or 
6 AU at 8.0 mm, and 58 mas or 7 AU at 9.8 mm. The errorbar 
at the bottom right indicates the systematic uncertainty on j^{R) 
resulting from absolute flux scale uncertainty. 


thus correspond to the slope of the line that intersects 
the points {x = log(A), y = \og{Y,xB^{Tx))/B^{T)}. 
By obtaining a large number of random samples of 
T,x X B^{Tx)/B^{T) we populate the PDF of A/3(i?) at 
each radius. We find the best-fit value of A/3 at a radius 
R from the peak of this PDF, and we derive confidence 
intervals from the region of the distribution that contains 
68.3%, 95.5%, and 99.7% of all samples at equal proba¬ 


bility (as described in Section for the uncertainties of 
the model parameters). These steps are performed for all 
radii sampled by the data, in order to obtain a best-fit 
value and inferred constraints on /?(i?) = /3c + 3^I3{R) 
for CYTau and DoAr 25. 

Figures [TT] and [I^ present the observational constraint 
on /3 as a function of radius, obtained for the circumstel- 
lar disks of CY Tau and DoAr 25. The values of (3 allowed 
by these multi-wavelength observations are significantly 
different from (and below) the ISM value of the dust 
opacity slope (/3ism ~ 1.7) for i? < 60 AU in the case of 
CYTau, and R < 70 AU for DoAr25. Furthermore, a 
gradient on j3 with radius is found in both disks. This 
gradient can only be consistent with a constant value 
of /? when the uncertainties in /3(i?) are extended to 7a. 
We note that the uncertainty in the absolute flux density 
scale will potentially introduce an additional systematic 
offset in the constrained values of /3(i?). But because 
a fractional uncertainty in the flux density scale affects 
all radii equally at a particular wavelength, the overall 
shape and significance of the deviation in P{R) from a 
constant value are not affected by uncertainties in the 
flux density scale. Given the 10-15% uncertainty in the 
absolute flux density scales for the three telescopes used, 
we estimate that the level of the systematic offset intro¬ 
duced to be ~ 0.16 for CYTau and ~ 0.08 for DoAr25, 
indicated by a ver tic al er rorbar in the bottom-right cor¬ 
ner of Figures [TT] and This systematic uncertainty is 

smallest for DoAr 25, because observations of this disk 
were obtained at wavelengths that are further separated 
than for CY Tau, reducing the uncertainty due to the 
increase in the lever arm used to infer (3. 


5.3. Radial variation of the dust opacity under the 
assumption of grain growth 

As discussed in the introduction, the inferred changes 
in the dust opacity for CY Tau and DoAr 25 may arise 
fro m changes in g rain size, composition, or grain geome¬ 
try. Draine ( 2006|) explored the frequency dependence of 
the dust opacity for different materials which may have 
enhanced emissivity at long wavelengths and thus can 
create low values of /3 with small grains. He concluded 
that changes in composition are unlikely to explain j3 < \ 
and advocates for changes in the grain-size distribution 
to explain the common finding of /3 1 in unresolved 

observations of protoplanetary disks (for references see 
§1). On the other hand, changes in grain por osity can 


produ c e low values of f3 (see, e.g.. Figure 11 of Kataoka 


et al. (2014)). However, the low filling factor needed 


to reach /3 < 1 requires quite large fluffy grains (e.g., 
larger than meter-size for filling factors of 10“^ — 10“^). 
Given the complexity of exploring changes in composi¬ 
tion, geometry, and size of dust grains simultaneously, 
and given that the option of large grains can more sim¬ 
ply reproduce the observed changes in kx, we assume 
that the radial variations of the dust opacity (observed 
as changes in /3 in the previous section) are only caused 
by changes in amax with radius. We now estimate the 
range of allowed maximum grain sizes at each location 
of the disk for CYTau and DoAr 25. 

As discussed by Perez et al. (2012), the determination 
of f3{R) depends on the assumption that the dust opacity 
follows a power-law with wavelength: kx oc . This es- 
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sential assumption might not be the case, as non-power- 
law dependence with wavelength in k\ can be seen on 
dust populations which are limited to a maximum grain 
size ( am.nT.) of sub- millimeter and mm-sized particles (see 
e.g., Draine|[2006 1. Thus, variations in the dust proper¬ 
ties of these disks, particularly Umax, should be directly 
constrained without this assumption. 


Following the method outlined by Perez et al. (2012), 
we find the best-fit dust opacity (which depends on 
dmax), along with the best-fit surface density, that is 
needed to reprod uce t he constraints on x x B^{Tx) 
found in Section |5.1[ Since the temperature profile at 
each wavelength [Tx) is found to be similar at all wave¬ 
lengths, we assume the dust temperature is equal to the 
mean Tlat e ach radius (see dashed line on left panel of 
Figures [9 and [l0|). At this point, we include the uncer¬ 
tainty in the absolute flux density scale in our modeling, 
which ranges between 10-15% for the three telescopes 
used. This improvement to our method is important, as 
it introduces another source of error in the constraints 
obtained from the flux density at each wavelength, and 
impacts the error budget in the characterization of ra¬ 
dial variations of the maximum grain size. We note that 
the uncertainty in the flux density scale as a function of 
wavelength wa s not included in o ur previous analysis of 
AS 209 (e.g., Perez et al. 2012). For each MGMG re¬ 
alization that reproduces the disk surface brightness as 
a function of radius, we draw a random number from a 
normal distribution centered at 1.0 and whose standard 
deviation (tr) corresponds to the flux scale uncertainty. 
Depending on the wavelength, this would correspond to 
tr = 0.1 or a — 0.15. We then scale the surface density 
normalization of each MCMC realization by a different 
random number, ensuring that the uncertainty in the 
absolute flux density scale is reflected in the resulting 
probability distribution. 

Assuming the same dust composition as in Sectionj^we 
place 3tT confidence limits for amax{R) and E(i?), while 
exploring two representative grain size distributions: a 
steep distribution with q — 3.5 and a shallow distribu¬ 
tion with q = 2.5. The chosen distributions are moti¬ 
vated by the resulting si ze-distribution of a collisional 
cascade ( ]Dohnanyi | 1969| and the measured grain size 
distribution in the ISM ( [Mathis et al. 1977), which are 
both characterized by n(a) oc while a less steep 

distribution of q = 2.5 is selected to observe the effects 
of a reduced total number of smaller particles. These 


results are presented as the shaded regions in Figures 13 
and 14 for CYTau and DoAr 25, respectively. We fim 
that within ~ 100 AU from the central protostar, grains 
have grown up to at least up to ~ 0.4 mm for GYTau 
and ^ 0.2 mm for DoAr 25. The best-ht amax{R) dis¬ 
plays a gradient with radius, where smaller grains are 
present in the outer disk and larger grains are present in 
the inner disk. However, the rise of Umax towards smaller 
radii strongly depends on the assumed slope of the grain 
size distribution. This occurs because a shallow grain- 
size distribution will have more large grains (by number) 
than a steep distribution. Since large grains dominate 
the emission at these wavelengths, we can reproduce the 
observed emission with a maximum grain size that can 
be smaller for a shallow distribution than for a steep dis¬ 
tribution. This explains the difference in our a^ax con¬ 


straints for the representative slopes q = 2.5 and q = 3.5 
that we selected. Unfortunately, current observations do 
not have the necessary sensitivity to discern the value of 
the grain size distribution slope. 

6. DISGUSSION 
6.1. Comparison with theory 


Birnstiel et al. (2012) presented simple analytical forms 
for the variation of Omax with radius for two different 
regimes: when a population of dust grains is limited by 
fragmentation of dust particles, whose differential veloc¬ 
ity is driven by turbulence in the disk: 


frag 


(i?) 


2 ^gas{R) 


^frag 


37r 


pat 


m 


( 6 ) 


and when a population of dust grains is limited by radial 
drift, as larger particles are removed from the population 
and accreted onto the star due to gas drag: 


dlnP 


d\nR 


( 7 ) 


drift / nt 2 l^dustjR) Vk{R) 

Here Jigaa and S^ust are the gas and dust surface den¬ 
sity profile, p is the dust grain internal densi ty, at is the 
disk turbulence parameter as presented by [Shakura fc| 
Sunyaev (1973), Cfrag is the velocity at which a collision 
between dust grains will result in fragmentation rather 
than growth, vk is the local Keplerian velocity, and Cs is 
the local sound speed of the gas. 

Employing our constraints on the disk surface den¬ 
sity, which assume a constant gas-to-dust ratio of 100:1 
over the entire disk, we compared our derived observa¬ 
tional constraints on Omax vs. radius with these grain 
growth models, adopting standard values for the turbu- 
lence (at = 0.01), du st density {p = 1.3 g cm“^, from the 
Pollack et al. (19941 abundances adopted in this paper) 


an d frag mentation velocity (cfrag = 10 m s 


aL][M2t . Also, we verify that mm-sized particles are m 


Birnstiel et 


the Epstein regime throughout the disk (with a Stokes 
number St < 1), which is a necessary condition for the 
use of equation!^ as the size limit in radial drift. For this 
particular set of turbulence and fragmentation parame¬ 
ters, the theoretical Omax curves for GY Tau and DoAr 25 


are shown in Figures 13 and |14[ respectively, where the 
continuous line illustrates the maximum grain size for a 
drift-dominated population, and the dashed blue line cor¬ 
respond to a fragmentation-dominated population. Note 
that since our constraint on the surface density with ra¬ 
dius depends on the adopted grain size distribution slope, 
q, the drift and fragmentation barriers will depend on q 
as well. 

In both disks and with this standard set of parameters, 
the fragmentation barrier under-predicts the maximum 
grain size achievable. Following equation!^ a better cor¬ 
respondence could be found if: (1) the disk turbulence 
were lowered to a* ~ 10“^ — 10“"*, or (2) the gas-to-dust 
ratio were higher than 100:1 by an order of magnitude, 
or (3) the critical velocity for fragmentation were higher 
than Cfrag = 10 m s“^. This last option seems unlikely 
if grains are compact, given that threshold velocities for 
fragmentation and/or bouncing of up to a few m s“^ 
have been measured experimentally for compact centime- 
ter and decimeter-sized grains (see e.g., Deckers & Teiser 
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Fig. 13.— Constraints for the CYTau disk on the surface density {top panels) and maximum grain size {bottom panels), for two grain 
size distributions q = 3.5 {left panels), q = 2.5 {right panels). The best-fit S(-R) (top) and amax{R) (bottom) are indicated by a continuous 
black line, while the shad ed region represents o ur confidence interval at 3cr. We compare our observational constraints with theoretical 
grain evolution models by|Birnstiel et al.|(|2012|| that include fragmentation (dashed lines, bottom panels), or radial drift (blue continuous 
line, bottom panels). For ttie tragmentation barrier we explore two different levels of turbulence in the CYTau disk: at = 0.01 (blue dashed 
line) and at = 0.001 (grey dashed line), with Uf^ag = 10 m s“^ in both cases. The smallest spatial resolution probed at each wavelength is 
indicated by a vertical dashed line in the top panels (at 7.1 mm, the smallest spatial scale corresponds to 42 mas or 6 AU, outside of the 
range plotted in this figure and indicated by an arrow). 


e for silica te grains, Heifiel mann et al. (|201 0 ) and 
al. (2015) for icy grains, and Blum fc Wurm |2008 1 


for a re^ew). On the ot her hand, it grains are porous 
(e.g. Kataoka et al. ]m4t , fragmentation could occur at 
higher threshold velocities. However, porous grains have 
larger cross sections resulting in mass loss through colli- 
sional erosion rather than fragm entation, which si milarly 
halts the growth of these grains (Krijt et al. Now, 

the second option of a higher gas-to-dust ratio lieeds to 
be pursued with measurements of the total gas mass in 
these disks in particular, but this alternative also seems 
unlikely given the recent results in the Taurus-Auriga 
star forming region where low gas-to-dust ratios (well 
below the ISM standard of 100: 1) are inferred for a larg e 
number of circumstellar disks by Williams & Best (2014). 
Finally, the first option of a quiescent disk of low tur¬ 
bulence {at < 0.01) for the disks around CY Tauand 
DoAr 25(and also AS 209) seems plausible, and would 
further augment grain growth. We test this option by cal¬ 
culating the expected maximum grain size when the disk 
turbulence is an order of magnitude lower {at = 0.001), 
and we hnd a better correspondence with our Umax con¬ 


straints on CY Tau and DoAr 25 (for this last disk we 
also increased the fragmentation velocity by 25%). Note 
that values of at ^ 10“^ — 10“^ would still produce tur¬ 
bulent relative velocities high enough for fragmentation, 
a condition attained when at > and confirmed 

for both CYTau and DoAr 25. The prediction of low 
turbulence should be tested with future observations of 
the turbulent linewidth in these circumstellar disks. 

At the same time, the radial drift barrier roughly 
agrees with our amax inference in both disks. For 
DoAr 25, the growth barrier imposed by a radial drift 
of macroscopic particles is within a factor of a few of the 
observational constraint, and the same occurs for CYTau 
at least up to 80 AU. Most likely it is both mechanisms 
that are at play simultaneously in these disks. Finally, we 
note that at a radii > 100 AU, the population of grains in 
CYTau is above both barriers by a factor of few at least, 
and this is the case for both a shallow or a steep grain 
size distributions. This result is readily appreciated in 
the CYTau data, where millimeter-wave emission at 1.3 
and 2.8 mm is observed at large distances from the star. 
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DoAr 25 (q=3.5) DoAr 25 (q=2.5) 
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Fig. 14.— Constraints for the DoAr 25 disk on the surface density {top panels) and maximum grain size {bottom panels), for two grain 
size distributions q = 3.5 {left panels), q = 2.5 {right panels). The best-fit S(i?) (top) and amax{R) (bottom) are indicated by a continuous 
black line, while the shad ed region represents o ur confidence interval at 3cr. We compare our observational constraints with theoretical 
grain evolution models by|Birnstiel et al.|(|2012f that include fragmentation (dashed lines, bottom panels), or drift (blue continuous line, 
bottom panels). For the fragmentation barrier we explore two different levels of turbulence and fragmentation velocity in the DoAr 25 disk: 
at = 0.01 with Ufrag = 10 m s“^ (blue dashed line), and at = 0.001 with = 12.5 m s“^ (groy dashed line). The spatial resolution at 
each wavelength is indicated by a vertical dashed line in the top panels (FWHM for 8.0 and 9.8 mm, HWHM for 2.8 and 0.9 mm). 


indicating the presence of mm-sized grains at ~ 100 AU. 
If the gas surface density distribution were different from 
the dust surface density we measured, the pressure gra¬ 
dient term | | in equation would change, and thus 

a region of higher gas pressure, i.e., a local maximum 
in P{R), could trap mm-sized dust grains in the outer 
disk of GYTau. In such a scenario, grains can overcome 
the radial drift barrier that is expected when the gas dis- 
tribution just monot onically increases with radius (e.g., 
Pinilla et al. 2012b). The origin of these pressure en- 
hancements is strongl y debated: locally they may arise 


from turbulent eddies (Klahr &: Henning|1997 Johansen 


& Klahr 20051, while on a global disk scale they may 


1999 


meria 


emerge from Kossby -wave instabilities ( Lovelace et al. 


Li et al. 20001, long-lived vortices (Large & Som- 
or from the gap opened by one or multiple 


planets (Pinilla et al. 2012a). However, higher angular 
resolutioh observations of the distribution of small and 
large grains, as well as observations of tracers in the gas 
emission, are needed to understand what is the likely 
cause of these large grains still present in the outer disk 
of GYTau. 


6.2. Comparison with other disks 

Over the past several years there have been multiple 
examples of circumstellar disks where the interpretation 
of radial variations of the dust properties is suppor ted by 
resolved observations at mm and cm wavelengths (Isella 
et al.||2010[ Guilloteau et al.||2011[ ITrotta et ah] 201J[ 


Menu et aL 20141 ^won et a h 2015p1 For example, iii 

'I'W Hya, |Menu et al. 'i2014() found that a radially ho- 
mogeneous distribution oFsinall and large grains (larger 
than ^ 100 pm) was not a suitable representation of the 
observations, and a steeper S(i?) for the larger grains 
was req uired to fit these mu lti-wavelength dataset. And 
recently |Kwon et al.| (|2015| studied a sample of disks 
and found that those with a steep midplane density gra¬ 
dient generally possess a smaller value of /3, indicative 
of dust growth and radial evol ution. The largest sampl e 
to date has been gathered by |Guilloteau et al.| (|201ip, 
where dual-wavelength observations of disks in 'iaurus- 
Auriga at angular resolutions between 0.4" — 1.0" were 
presented. This study found evidence for changes in the 
value of /3 for several sources including CY Tau. Albeit 
limited by the angular resolution of the lowest frequency 
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Fig. 15.— Left: Observational constraints on the spectral index of the dust opacity /3, with Ic r error bars, as measured from published 
resolved millimete r observations, where all but CQ Tau have been observed only at 1 and 3 mm (|Isella et al.|2010| |Guilloteau et al.|2011| 
|Trotta et al.|2013f. Note that the majority of the previously available constraints are still consistent with a constant value ot p across the 
disk at '6a. Ihe rnean profile of fi{R) based on these observations is shown by the gray shaded region. Right: The radial variat ion of /3 for 
the three circumstellar disks analyzed in a consistent fashion: CYTau and DoAr25 (this paper), and AS 209 (|Perez et al.|2012||, where the 
shaded regions show the la- constraints. We note that at the 3 (t level the 3 disks analyzed here at consistent with each other and with the 
average /3(R) profile from the literature. 



band, which corresponds to ^ 100 AU resolution for their 
2.7 mm observations of CY Tau, these authors infer a 
change from /3 ~ 0 in the inner disk to high /3-values 
in the outer disk of this young star. However, given the 
closeness in frequency o f the two wavelengths studied by 
Guilloteau et al. (2011), together with the lack of high 
spatial resolutioii at 2.7 mm, implies that their constraint 
is consistent with a constant value of f3 at th e 3a level for 
this particular object. As shown in Section 5.21 our in¬ 


creased wavelength coverage and high sensitivity allows 
us to discriminate radial variations of the dust emissivity 
index for CY Tau with high significance, where a con¬ 
stant value of /3 is excluded at more than 7a. This arises 
because the uncertainty in the spectral index {aa) is di¬ 
rectly related to the wavelength coverage and SNR of 
the observations. In the case of dual-wavelength observa¬ 
tions, simple error propagation results in an uncertainty 
^2 _ 


ofcr„ = (ln(i/i/i/ 2 )) ^(SNR,^ -I- SNR ), which indicates 


i-2^snr-2) 

that an efficient way to reduce the uncertainty in a, and 
thus in /3, is to expand the wavelength coverage. 

To compare radial changes in the dust properties for 
different circumstellar disks, we have gathered measure¬ 
ments from the literature that attempt to constrain /3(i?), 
mostly from observations at the 1 and 3 mm bands. Since 
most of the available measurements barely resolve each 
circumstellar disk, we employ only two representative 
measurements of /3, one from the inner disk regions and 
another from the outer disk. Thus, for each published re¬ 
sult we compile the value of /3 and its Icr uncertainty at 
the best resolution of the observations (between ~ 20—50 
AU depending on each target) to illustrate the inner disk 
constraint in /3, and at either the observed outer disk 
radius or at i? = 120 AU, whichever is smallest, to rep¬ 
resent the constraint of 13 in the outer disk (except for 
RY Tau and DC Tau, which at the outer disk radius have 
an unconstrained value of /3; for these we select the out¬ 
ermost radii for which a /3 constraint is still significant). 

The compiled literature measurements are presented 
in the left panel of Figure [TS} where we have computed 
the slope in /3(i?) between the inner and outer disk for 


each star (dashed lines), and the average slope in j3{R) 
for all the disks analyzed in the literature (grey shaded 
region). The steepest change in the dust emissivity spec¬ 
tral index occurs for DC Tau B (where /3 = 0.1 ± 0.3 at 
20AU, incr easing to j3 = 1.5 ± 0.3 at 120 AU; |Guilloteau| 
et al.|20lT ). Excluding the added uncertainty in /3 due to 
the absolute flux scale uncertainty for the two millimeter- 
wave bands observed for DG Tau B, the resulting change 
in P{R) for this object is comparable to the constraints 
found for GY Tau, DoAr25, and AS 209, which are pre¬ 
sented together on the left panel of Figure [TS] 
Interestingly, the mean /3(i?) profile from literature 
measurements is less steep than for the 3 disks analyzed 
in a consistent fashion. This could be due to the fact 
that most previously published observations only include 
dust emission traced at millimeter wavelengths (1 and 
3 mm bands), which may be optically thick and thus 
reduce the overall value of /3. Since the young stars pre¬ 
sented in the left panel of Figure boast some of the 
brightest disks observed in millimeter continuum emis¬ 
sion, they represent a population of massive disks where 
high optical depth at millimeter wavelengths may be ex¬ 
pected. Thus, including long-wavelength observations is 
critical for the calculation of P{R), not only for the in¬ 
creased lever arm but also to avoid regions of high op¬ 
tical depth. We note that the comparison between pre¬ 
vious results in the literature (/3nim(3?)) and our mul¬ 
tiwavelength analysis (/3mm — cm(i?)) should ideally be 
performed over a similar wavelength range to avoid po¬ 
tential biases between different wavelengths that trace 
different grain sizes, and could explain some of the dif¬ 
ferences between /3mm (3?) and /3mm-cm(3?) seen in Figure 
[Tsl Finally, we note that GY Tau shows the steepest rise 
in /3 with radius, indicative of a substantial change in 
dust properties over tens of AU in this disk, while the 
smoothest rise in the emissivity spectral index occurs for 
AS 209 in the inner disk and for DoAr 25 in the outer 
disk. However, these differences at not significant when 
the disks are compared at the 3tT level. 
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7. CONGLUSIONS 

We present multi-wavelength observations of the proto¬ 
planetary disks surrounding the young stars GY Tan and 
DoAr 25. Observations at 0.9, 1.3, 2.8, 8.0 and 9.8 mm 
spatially resolve the dust continuum emission from these 
disks, down to scales of tens of AU. From observations 
obtained at 5 cm, we quantify the amount of continuum 
emission whose origin is not thermal dust emission. We 
estimate the level of contamination for DoAr 25 to be 
< 24% of the integrated emission at 8.0 and 9.8 mm, 
while for GY Tau we establish an upper limit of con¬ 
tamination of < 4% of the total emission at 7.1 mm. 
Although these levels of contamination may seem small, 
they need to be quantified when characterizing the dust 
emission spectrum. The temperature and surface density 
profiles for both disks were constrained at all observed 
wavelengths shorter than 5 cm, and from this modeling, 
we hnd that a constant dust opacity with radius does 
not fit these observations. This result is supported by 
the observational fact that the normalized visibility pro¬ 
files of the dust continuum emission differ for different 
wavelengths: cm-wave observations trace emission from 
a compact disk structure, while mm-wave observations 
trace a more extended disk structure. 

Since our modeling indicates that the dust continuum 
emission is optically thin over the spatial scales explored 
(and this is also warranted by the observed brightness 
temperature prohles in §3.4), the observed emission spec¬ 
trum will be directly proportional to the dust opacity 
spectrum, which allows us to infer radial variations of 
the dust opacity spectral index, We find that a 

constant value of /? is excluded by our observations at 
more than 7 ct in both disks. Close to the central pro¬ 
tostar, for i? < 50 AU, the constraints on /3(i?) are of 
high significance and indicate a rapid change of the dust 
properties between the inner disk and the outer disk. In 
the outermost regions of the disk where the dust emis¬ 
sion has lower SNR, we can only set an lower limit for 
the dust opacity spectral index and we find that /3 > 1.0 
for CYTau at ii > 80 AU, while /3 > 1.2 for DoAr 25 at 
R> 80 AU. We find that for the disks of AS 209, CY Tau, 
and DoAr 25, the profiles of f3{R) are steeper than for the 


level dlsella et al. 

2010 

Guilloteau et al. 1120 11 

Trotta et 

al. 

2U13), while t 

hese are all consistent at the '6a level. 


by 

of 


Perez et al. (2012) is critical for obtaining constraints 
ligh significance on (5{R). 


Assuming that the changes in the dust properties dis¬ 
cussed above arise only from changes in the maximum 
particle size of the disk grain size distribution, we derive 
radial variations of the maximum grain size, ttmaxiR)- 
For the two different grain size distribution explored, 


q = 2.5 and q = 3.5, we find a gradient in Umax with 
increasing grain size at smaller radial distance from the 
star, in both the CY Tau and DoAr 25 disks. We compare 
the observational constraints in amax{R) with theoretical 
expectations of the maximum grain size when a popula¬ 
tion of grains grows and evolves while being limited by 
fragmentation or by radial drift. For an assumed disk 
turbulence of at = 0.01, the fragmentation barrier does 
not seem to be the limiting factor in the growth of these 
grains, as we observe millimeter and centimeter dust par¬ 
ticles at distances from the star where the fragmentation 
barrier should have depleted these large particles. How¬ 
ever, our observations seem to be consistent with a low- 
turbulence disk midplane with at ~ 0.001 in CY Ta u 
and DoAr 25, as well as in AS 209 (Perez et al. 2012). 
If disks are generally low-turbulence, then fragmentation 
may play a lesser role in shaping the grain size distribu¬ 
tion of protoplanetary disks. Future constraints on disk 
turbulence will be necessary to assess this process. For 
the CY Tau disk, the population of grains inferred from 
our observations is above the growth barrier imposed by 
radial drift and fragmentation, with a difference of of a 
factor of a few at > 100 AU radius. We hypothesize 
that this difference could be readily explained if a re¬ 
gion of higher gas pressure exists in the outer disk, and 
is keeping these millimeter-sized particles from drifting 
inwards. Future high resolution observations that trace 
millimeter dust grains, as well as the gas distribution, 
are needed. Observations with the Atacama Large sub- 
Millimeter Array (ALMA) will be extremely well suited 
for such investigations, as spatial resolutions matching 
the VLA observa tions at 7 mm (of 0.07" or bette r) are 
already possible (ALMA Partnership et al. |2015 ). Fi¬ 
nally, future observations of multiple disks at different 
stages of evolution (classical vs. transitional), as well as 
different stellar mass or disk masses, will provide valu¬ 
able information related to the first steps toward planet 
formation. 
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